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Abstract 

We present a new model describing strongly correlated electrons on a general 
ci-dimensional lattice. It differs from the Hubbard model by interactions of nearest 
neighbours, and it contains the t-J model as a special case. 

The model naturally describes local electron pairs, which can move coherently at 
arbitrary momentum. By using an r^-pairing mechanism we can construct eigenstates 
of the hamiltonian with off-diagonal-long-range-order (ODLRO). These might help 
to relate the model to high-T^ superconductivity. 

On a one-dimensional lattice, the model is exactly solvable by Bethe Ansatz. 
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The study of strongly correlated elctrons on a lattice is an important tool in 
theoretical condensed matter physics in general, and in the study of high-T^ super- 
conductivity in particular. Two well-studied models are the Hubbard model and 
the t-J model. On a one- dimensional lattice these models are both exactly solvable 
by Bethe Ansatz. In this letter we propose a new model, which is again solvable in 
one dimension, and which combines and extends some of the interesting features of 
the Hubbard model and the t-J model. 

Electrons on a lattice are described by operators Cj^a , j = ^, ■ ■ ■ , L, a =t, i, 
where L is the total number of lattice sites. These are canonical Fermi operators 
with anti-commutation relations given by {c|^,Cj^T-} = 5ij5o-,r- The state |0) (the 
Fock vacuum) satisfies Ci^„\ ) = 0. At a given lattice site i there are four possible 
electronic states 

|0), IT), = 4^10), |i), = cljO), \nh = clcl\0). (1) 

By ^j o- = cJ^Cj^o- we denote the number operator for electrons with spin a on site i 
and we write = nj ^ + n^^j^. The spin-operators S = Ylj=i Sj , , S'^ , S^, 

Sj = cl^Cj^i , S] = 4^Cj- 1 , S] = ^(n^- 1 - n^- 1) , (2) 

form an SU{2) algebra and they commute with the hamiltonians that we consider 
below. [We shall always give local expressions Oj for symmetry generators, implying 
that the global ones are obtained as O = Y.j=i Oj.] 
The Hubbard model hamiltonian can be written as 

^Hubbard = _Y.T. {^Wo + etc,,.) + t/E(^.,T - - \) > (3) 

where the first summation runs over all nearest neighbour pairs (jA;). It contains 
kinetic (hopping) terms for the electrons and an on-site interaction term for electron 
pairs. An interesting feature (on a bipartite periodic lattice) is the so-called ?7-pairing 
symmetry [||, H, which involves operators and rfn which form another SU{2) 
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algebra, and which commute with the hamiltonian (j^). Using this symmetry one can, 
starting from an eigenstate of the hamiltonian, create a new eigenstate r/J^I?/'), 
which contains an additional local electron pair of momentum vr. The spin SU{2) 
algebra @ and the ?7-pairing SU{2) algebras together form an 5*0(4) symmetry 
algebra. In one dimension, the Hubbard model is solvable by Bethe Ansatz [Q. 

In the t-J model, there is a kinematical constraint which forbids the occurence 
of two electrons on the same lattice site. On this restricted Hilbert space the t-J 
hamiltonian (with t = 1, J = 2) acts as = — J2{jk) Hj~k^^ with 

-2 (5|5,^ + ^{S]S, + S,Si) -\{l-n,- n,) - ^n,n,^ , (4) 

where we defined 

<3i,T = (1 - c^,^ , Qix = (1 - nj^{) Cj^^ (5) 

and the operators Sj, Sj and Sj are as in (0). The t-J model (§) is supersymmetric 
and the spin SU{2) symmetry algebra gets enlarged to the superalgebra SU {1\2) 
0, (see for the description and classification of the classical Lie superalgebras) . 
The generators of this symmetry algebra are S", S\ S^, Q-^, Qp Qi, Q\ and T = 
2L — J2j=i ^j. In one dimension the supersymmetric t-J model (|) is exactly solvable 
by Bethe Ansatz [0, g |. 

Before we present the hamiltonian of the new model, we give some motivation, 
which is based on what we know about the materials that exhibit high-T^ super- 
conductivity. It has been found that the electrons in these materials form 'Cooper 
pairs', which are spin-singlets, and that these pairs are much smaller than in the 
traditional superconductors. As a limiting case one can consider models which have 
electron pairs of size zero, i.e., pairs that are localised on single lattice sites. We 
will call such localised electron pairs localons. 

In the t-J model localons are ruled out by the kinematical constraint on the space 
of states, and in the Hubbard model only local pairs of momentum vr exist. Below 
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we shall see that in our new model localons can move coherently with arbitrary 
momentum. Apart from these local pairs, the new model may also have bound 
states that are finite-size electron pairs. 



Let us now present the hamiltonian of the new model on a general ci- dimensional 
lattice. We write it as 



H = H' + U - ^)inu - 3) + E + ^ ' . (6) 

j=i j=i j=i 



where is given by 



H° = -T.Hlk, aren.n. (7) 



with 

Hlk = 4,TCi,T(l - "-i.i - "fc,i) + cIt^^.tII - - ^kd) 

- 1)K - 1) + clAd^kdCk,^ + Ci,iCi,T4,T4a 

11 11 

+ Kt - 2^if^j,l ~2^ + - 2^irik,i - -) . (8) 

This hamiltonian contains kinetic terms and interaction terms that combine those 
of the Hubbard and of the t-J model. The second term in is the on-site Hubbard 
interaction term (notice that it also gets a contribution from H^). The third and 
the fourth terms in (|^) introduce a chemical potential fi and a magnetic field h. 
Roughly speaking, the new model can be viewed as a modified Hubbard model with 
additional nearest neighbour interactions similar to those in the t-J model. 

The hamiltonian is invariant under spin refiection Cj | ^ Cj | and under 
particle-hole replacement c] ^. ^ Cj^a- In addition to the spin SU{2) generators (|]), 
the following operators commute with 
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?7-pairing SU{2), with generators r], ri\ r]^ 



1 1 



Together with the spin SU{2) algebra this gives an 50(4) algebra which 
is similar to the one for the Hubbard model. This symmetry makes it possible 
to generalise the r^-pairing mechanism, which was developed for the Hubbard 
model in p|, to the new model. 

• Supersymmetries. There are eight supersymmetries in total: Q^, Q|, Qi and 
Q| given in (H) and the operators Qo-; Qt^ 

QjA = %,iCj,T . Qia = %,TCj,i ■ (10) 

These generators, together with the operator SjLi 1 (which is constant and equal 
to L), form the superalgebra SU {2\2). [Like S'f/(4), this algebra has 15 generators, 
eight of which are fermionic. In the fundamental representation, the generators can 
be represented as 4 x 4 supermatrices with vanishing supertrace 0.] 

The symmetries of the hamiltonian can be made manifest as follows. We 
first add one more generator to the symmetry algebra, which is 

X = j:X,, X, = (n,-t-b(%i-^) , [H\X] = Q. (11) 
i=i ^ ^ 

This extends the superalgebra SU {2\2) to U{2\2). We denote the generators of 
this algebra by J^, where a = 1,2, . . . , 16. We now introduce an invariant, non- 
degenerate 2- index tensor, denoted by K"'^, which is the inverse of Kais = str^JaJp), 
where the Ja are 4x4 supermatrices in the fundamental representation. Using this, 
we can cast H^j^ in a group-theoretical form as follows 

16 

Hi = E K'^''J3,^Jk,P (12) 

a,(3=l 

~ E (Q],rTQk,(7 + Qk,aQj,<r ~ Q],aQk,cr — Ql:,crQj,<r) 

+ (2 r^^r^l + ri% + vA) " (2 S^^Sl + ^j^^ + S,Sl) + X, + . (13) 
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It is easily checked that this expression agrees with the formula for Hjf, in (|^). The 
expression ([T^) immediately makes it clear that if" commutes with all 16 generators 
of f/(2|2). 

We would like to stress that the appearance of the algebra f/(2|2) in the model 
is not too surprising: on each lattice site there are four electronic states (|l|), two 
of which are fermionic. The supergroup U{2\2) is the group of all unitary rotations 
of these four states into one another. Our hamiltonian has been chosen such 
that it commutes with the entire algebra U{2\2) and is therefore very natural. The 
analogous construction for U{1\2) leads to the t-J hamiltonian and for U(2) it 
leads to the spin-| XXX Heisenberg model. 

The spectrum of the hamiltonian is symmetric around zero. This follows 
from the discrete symmetry cjj^ ^ Cj^i, for which ^ —H^. 

There is a further aspect of that deserves to be mentioned: the terms Hj^. 
act as a graded permutations of the electron states (|^) at sites j and k. By 'graded' 
we mean that there is an extra minus sign if the two states that are permuted are 
both (fermionic) single electron states. For example, 

HlcU0) = cl^\0), ii°,4cI,JO) = -44^^10), etc. (14) 

In this respect, the new model generalises the spin-| XXX model and the t-J 
model (^). The nearest neighbour hamiltonians of these models have a similar 
interpretation as graded permutations of the basic states, which are {|T )? |i )} for 
the spin-| XXX model and {| ), || ), || )} for the t-J model. Lattice hamiltonians 
that act like (graded) permutations were first considered by Sutherland in p. 

We define the number operators N^, Nj^ (the number of single electrons with 
given spin) and A*"; (the number of localons) by 

L L L 

N^ + Ni = Y^ ^i,T . N^ + Ni = Y^ n,, i , Ni = J2 ^.,T (15) 
j=i j=i j=i 

and we write = + Nj^ for the total number of single electrons. The fact that 
is a permutation makes it clear that these number operators commute with H^, 
so that can be diagonalised within a sector with given numbers A'^^, and Ni. 
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This implies that the terms proportional to U, fi and h in (^), which break the 
symmetry f/(2|2), will not affect the solvability of the model in one dimension. 

In the sectors without localons reduces to the t-J hamiltonian (|^). (This is 
clear from the fact that they both act as permutations.) The new model reduces to 
the spin- 1 XXX model in the sector with only vacancies and localons, and similarly 
in the (half-filled) sector with one single electron on each site. 

Let us now briefly discuss some physical aspects of the new model. We first 
remark that we can always (for general lattices in an arbitrary number of dimensions) 
construct a number of exact eigenstates of the hamiltonian which show off-diagonal- 
long-range-order (ODLRO), which is characteristic for superconductivity |Ty]. For 
this we follow the construction which was developed for the Hubbard model by C.N. 
Yang in [§. The state = (^^)^|0) is an eigenstate of the hamiltonian with 
energy E = + — M, where M is the total number of nearest neigbour 
links {jk) in the lattice. Following f^, we compute the following off-diagonal matrix 
element (A; 7^ /) of the reduced density matrix p2 

{{k, i){k, t)| P2 1(^ T)(/, i)) = , , , = T77 — TV • 

The fact that this off-diagonal matrix element is constant for large distances \ j — k\ 
establishes the property of ODLRO for the states \I/Ar. In the thermodynamic limit 
we have ODLRO as soon as the occupation number N of the zero-momentum localon 
state becomes macroscopic, i.e., when N/L —>■ D becomes finite. 

The local electron pairs that participate in the r^-pairing have momentum zero. 
However, the model also admits localons that move with arbitrary momentum. 
This follows from the fact that Hj;^ acts as a permutation of the electronic states 
(|I]) on neighbouring sites: due to this localons cannot decay and move coher- 
ently. On a c?- dimensional square lattice (with lattice spacing a) the wavefunction 
J2x exp(?a;fc) ct^ct J 0), which describes a single localon of momentum k over the 
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bare vacuum, is an exact eigenfunction of the hamiltonian (^) of energy 

d ^ 

E = 2d-2Y, cos{kma) + f/- + 2/x - M . (17) 

m=l ^ 

Multi-localon wave-functions, as well as wavefunctions with single electrons, exist 
but cannot easily be written down for higher-dimensional lattices. However, in 
one dimension the model is exactly solvable by Bethe Ansatz (BA), and we can 
obtain explicit expressions for general eigenstates of the hamiltonian. We think that 
it is worthwhile to study this exact solution, and that this will lead to a better 
understanding of the higher-dimensional model as well. 

We will here briefly summarise the results of the exact solution in one dimension; 



the details are deferred to a separate publication ||TT|. The exact solution starts from 
the observation that the hamiltonian is a graded permutation ([l^ of four states, 
of which two are fermionic and two are bosonic. The BA analysis for hamiltonians 
which are graded permutations was first considered by Sutherland in P, see also 
|T^. The method of solution is the algebraic version of the 'nested Bethe Ansatz' 



|[T3[| (for an introduction to the algebraic BA, see [|T4|)- Each step of the nesting 
involves the introduction of a set of spectral parameters, which are in our case Xj, 
Xl^^ and Ap \ where j = 1, {N^ + Ni), k = 1, and I = 1, ... , N^. 

For each choice of a set of rapidities we can construct an eigenstate of the hamil- 
tonian in the sector specified by Ni, Ne and Ni, with energy given by 

E'= E XT-T-L. (18) 

j=l + 4 

The boundary conditions for these eigenstates lead to the following set of Bethe 
equations for the rapidities Aj, A^^"* and Ap"* 

VA. + ty M + ^ iiiA«-A, + f 

n ^k - ^j' + 2 _ TT - '^U +2 

xW X i x(i) x(2) i 

j'=l ^k - ^i' ~ 2 ''=1 ^k ~ \' ^2 
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n 



1 ^i' ~ ^, 



(2) 



n 



1 '^k' 



X 



(2) 



(19) 



These equations, together with the expression ( [I8| ) for the energy, guarantee that 
we shall be able to describe explicitly the ground state of our model for arbitrary 



density of electons, coupling constant U and magnetic field h [O . 
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